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根 肿 菌 侵 染 下 窗 蓝 生物 碱 合成 机 制 
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摘 要 : 为 了 探究 根 肿 菌 胁 迫 对 钦 蓝 生物 碱 及 其 合成 关键 酶 基因 表达 的 影响 , 该 文 对 根 肿 菌 
侵 染 后 0、7、14、21 d 的 获 蓝 进行 病情 形态 分 级 、 组 织 学 观察 、 生 理 生 化 指标 测定 以 及 转 
录 组 学 和 代谢 组 学 分 析 。 结 果 表 明 : (1) 接 菌 后 0、7、14、21 d 杰 蓝 根部 分 别 发 展 为 0 级 、 
1 级 、3 级 、5 级 的 肿 根 ， 且 7d 是 皮层 入 侵 的 关键 时 间 点 。(2) 接种 根 肿 菌 14 d Ja, PH 
叶 内 可 溶性 蛋白 、 丙 二 醛 含 量 、 超 氧化 物 歧化 酶 、 过 氧化 物 酶 、 多 酚 氧 化 酶 和 过 氧化 氢 酶 活 
性 与 同时 间 对 照 组 比较 显著 提高 ， 并 随 着 接 菌 时 间 的 延长 旦 增加 的 趋势 。 (3) 代谢 组 学 一 
检测 到 161 种 生物 碱 ， 其 中 呀 唆 类 生物 碱 数量 较 多 ; 与 未 接 菌 相 比 ， 获 蓝 接 菌 后 7、14、 
21 d 分 别 存 在 16、17、39 种 差异 代谢 物 且 各 组 差异 代谢 物 多 富 集 在 生物 碱 和 氨基 酸 代 谢 通 
Ko (4) 转录 组 测序 结果 显示 : 与 未 接 菌 相 比 ， 徐 蓝 接 菌 后 7、14、21 d 分 别 存在 2 439、 
256、6 437 个 差异 表达 基因 ， 这 3 组 共同 富 集 到 11 个 生物 碱 相 关 的 代谢 通路 ， 与 未 接 菌 相 
比 ， 接 菌 后 7、14、21d4 有 9 个 基因 (编码 4 FHM THS, TAT, YUCCA, ALDH) 表达 量 均 
上 升 。 以 上 结果 揭示 了 芸 草根 肿 菌 与 黎 蓝 之 间 的 互 作 机 制 ,探究 了 根 肿 菌 对 叫 唆 生物 碱 合成 
及 其 关键 酶 基因 的 影响 ， 为 后 期 研究 徐 蓝 根 肿 病 抗 性 基因 及 生物 碱 次 生 代 谢 途 径 黄 定 基础 。 
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Abstract: To explore expression level of alkaloid and its synthetase gene in /satis indigotica upon 
Plasmodiophora brassicae Woronin exposure. Our study verified grade of disease severity 
according to morphology. Moreover, histological, physiological, biochemical parameters have 
been collected together with transcriptomics and metabolomic analysis in Isatis indigotica after 
infection at time scale 0, 7, 14, 21 d. Results showed that: (1) After infect with P brassicae, I. 
indigotica showed club roots grade in 0, 1, 3, and 5 respectively at 7, 14, 21 d, notably, cortical 
invasion occurred on 7 d. (2) Our study revealed when P. brassicae exposes lasting 14 d later, the 
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level of soluble protein and malondialdehyde, along with superoxide dismutase, peroxidase, 
polyphenol oxidase and catalase activity in Z. indigotica increased significantly compared to 
control group at time depended manner. (3) A total of 161 alkaloids were detected in 
metabolomics, among those alkaloids, indoles were noticed as the most abundant form. There 
were 16, 17 and 39 discriminating metabolites have been spotted after infected with P. brassicae 
at 7, 14, 21 d, the most discriminating metabolites enriched at alkaloid and amino acid metabolism 
pathways. (4) Transcriptome analysis showed that there were 2 439, 256 and 6 437 genes 
expression alteration at 7, 14, 21 d compared to control, those differentially expressed genes 
enriched at 11 alkaloids related metabolism pathways. Markedly, expression level of 9 genes 
(encoding for enzymes thebaine synthase, tyrosine aminotransferase, indole-3-pyruvate 
monooxygenase and aldehyde dehydrogenase) were increased after infection at 7, 14, 21 d. Our 
results revealed the interaction between P brassicae and I. indigotica. The effects of P. brassicae 
on indole alkaloids and enzymic genes expression laid ground for studying genes mediate club 
root resistance and alkaloid secondary metabolic pathways in Z. indigotica. 
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PAW Qsatis indigotica) ， 十 字 花 科 获 蓝 属 的 二 年 生 植 物 ， 全 国 各 地 均 有 栽培 ， 具 有 很 
大 的 经 济 和 药 用 价值 ， 其 干燥 根 、 干 燥 叶 为 清热 解毒 的 代表 药 板 蓝 根 和 大 青 叶 (Wong et al., 
2022) 。 窗 蓝 中 的 主要 生物 活性 成 分 生物 碱 具有 抗 病 毒 、 抗 菌 及 免疫 调节 的 药理 作用 ， 其 中 
叫 唆 类 生物 碱 的 药理 活性 报道 最 多 ， 已 经 证 实 靛 蓝 、 靛 玉 红 、 表 告 依 春 、5$- 羟 基 呀 唆 、 叫 吹 
-3- 甲 醛 和 色 胺 酮 具有 抗 病毒 、 抗 肿瘤 、 白 细胞 抑制 及 抑 菌 活性 (Sinha et al., 2008; Chen et al., 
2021; 杨 立 国 等 ，2021; ) 。《 中 国药 典 》 也 规定 板蓝根 和 大 青 叶 药材 及 其 制剂 质量 控制 的 
重要 指标 分 别 是 表 告 依 春 和 靛 玉 红 中国 药典，2020) 。 
十 字 花 科 根 肿 病 俗称 大 根 病 ， 其 病原 为 丢 草 根 肿 菌 CPlasmodiophora brassicae) ， 是 一 
种 土 传 病害 ， 专 门 寄 生 在 检 蓝 等 十 字 花 科 作 物 〈 和 白菜、 油菜、 葛 卜 、 芥 菜 ) 的 根部 ， 通 常 
侵入 根 毛 ， 再 产生 次 级 游 动 息 子 入 侵 皮 层 组 织 ， 致 使 皮层 细胞 增 大 ， 挤 压 变 形 ， 侵 染 后 期 细 
HA er AE RE, 在 此 过 程 中 植株 根部 肿 大 呈 瘤 状 且 生长 发 育 迟 绥 、 凋 萎 下 重 。 根 肿 病 
侵 染 也 会 诱导 寄主 抗 性 响应 如 产生 大 量 的 氧 自 由 基 ， 促 使 膜 脂 过 氧化 ， 破 坏 细胞 膜 的 结构 ， 
MIRAN (Catalase, CAT) 、 多 酚 氧 化 酶 (Polyphenol oxidase, PPO) 、 超 氧化 物 靶 化 
li (Superoxide dismutase, SOD) 和 过 氧化 物 酶 (Peroxidase, POD) 等 防御 酶 主要 参与 活 
性 氧 的 消除 ， 维 持 细胞 膜 稳 定性 〈Qi et al., 2020; Huang et al., 2022; EAI, 2022). Æ 
六 月 (2021) 研究 表明 抗 病 大 白菜 DH40R 中 的 CAT. SOD, POD 活性 和 可 溶性 蛋白 含量 在 
根 肿 菌 侵 染 8 d 明显 升 高 ， 强 烈 抵 抗 根 肿 菌 侵 染 ;朱红 芳 等 (2015) 研究 表明 : 初期 筛选 搞 
肿 病 的 品种 时 ， 将 可 溶性 糖 、 可 溶性 蛋白 和 MDA 作为 主要 指标 ， 郭 珍 (2018) 研究 发 现 
CAT 与 油菜 抗 病 性 相关 ，SOD、POD 与 油菜 抗 病 性 呈正 相关 。 
目前 , 对 十 字 花 科 植 物 感染 根 肿 菌 的 相关 研究 主要 还 是 集中 在 靶 荤 属 植 物 , 对 拟 南 芥 的 
研究 较 少 , 徐 蓝 感染 根 肿 菌 的 防御 机 制 以 及 体内 代谢 变化 的 研究 几乎 没有 报道 。 因 此 迫切 需 
要 探究 获 蓝 与 根 肿 菌 之 间 的 互 作 机 制 ， 以 及 根 肿 菌 对 获 蓝 有 效 成 分 积累 的 影响 。 本 研究 以 接 
菌 后 0、7、14、21 d 的 获 蓝 为 研究 对 象 ， 采 用 Humina 高 通 量 测序 技术 和 超 高 效 液 相 色谱 - 
串联 质谱 技术 ， 通 过 病情 形态 分 级 、 组 织 学 观察 、 生 理 生 化 指标 测定 以 及 转录 代谢 分 析 ， 拟 
探讨 以 下 问题 : C) 根 肿 菌 侵 染 获 蓝 的 动态 过 程 ，(2) 徐 蓝 抗 根 肿 病 的 生理 防御 机 制 ; (3 ) 
根 肿 菌 胁 迫 对 获 蓝 生物 碱 类 化 合 物 合成 机 制 的 影响 。 本 研究 能 够 为 深入 挖掘 获 蓝 抗 根 肿 病 功 
能 基因 及 次 生 代谢 物 合成 途径 等 分 子 研究 芮 定 基 础 。 


1 材料 与 方法 


1.1 材料 

供 试 菌株 和 黎 蓝 品种 : 菌株 采 自 昆明 市 盘 龙 区 阿 子 营地 区 的 大 白菜 发 病根 瘤 , 经 云南 省 
农业 科学 院 鉴 定 为 4 号 生理 小 种 (Williams 鉴定 系统 ) 。 实 验 所 需 小 叶 获 蓝 种 子 也 由 云南 省 
农业 科学 院 园 艺 作物 研究 所 的 十 字 花 科 课 题 组 张 丽 琴 老师 提供 并 鉴定 。 

试剂 及 仪器 : 台 盼 蓝 溶 液 (SIGMA ) 、 甲 醛 -醋酸 -乙醇 固定 液 〈50%) ; 牛 血清 蛋白 、 
核 黄 素 、 愈 创 木 酚 、2- 硫 代 巴 比 妥 酸 ， 北 京 索 莱 宝 科技 有 限 公 司 ; 邻 葵 二 酚 、 三 氧 乙 酸 溶液 、 
L- 甲 硫 氮 酸 ， 上 海 麦 克 林 生 化 科技 有 限 公 司 ; 氧化 硝 基 四 氮 唑 蓝 (NBT) 、 乙 二 胺 四 乙酸 
EDTA， 德 国 Biofroxx 生物 试剂 ， 过 氧化 毛 30%， 云 南 景 锐 科技 有 限 公 司 。 其 他 仪器 和 耗材 
包括 : Nikon 数码 显微镜 ， 血 球 计数 板 ， 烧 标 ， 离 心 管 ， 容 量 瓶 ， 量 简 ， 育 苗 盘 ， 研 钵 ,而 
棒 ， 比 色 下 ， 普 析 TU-1810 紫外 可 见 分 光 光 度 计 ， 超 低温 冷冻 离心 机 ， 光 照 培 养 箱 ， 台 式 
喜 风 干燥 箱 ， 电 热 恒温 水 浴 锅 ， 千 分 之 一 电子 天 平 ， 托 盘 天 平 ， 冷 冻 冰 箱 等 。 
1.2 试验 处 理 

将 冷冻 的 大 白菜 发 病根 痛 置 于 室温 解冻 ， 采用 杨 佩 文 等 (2002) 的 方法 提取 根 肿 菌 休眠 
FF, 血球 计数 板 测 定 根 有 种 菌 孢子 甚 浮 液 的 浓度 , 稀释 成 5x107 CFU mL! 病原 菌 菌 液 备 用 。 
JH RSE 72 孔 育 苗 盘 中 栽培 ，1 个 穴播 种 3 颗 种 子 ， 待 其 长 出 2~ 3 对 真 叶 后 间苗 ， 注 射 
法 进行 根 肿 菌 (1x10” CFU mL! 休眠 孢子 ) 接种 ， 用 10 mL 注射 器 吸取 菌 液 ， 将 菌 液 缓慢 
注射 到 实验 苗 根部 ,2 ml 菌 液 / 株 。 分 别 在 接种 根 肿 菌 后 0、7、14、21 d 采 和 集 正 常生 长 (BLG-CK1、 
BLG-CK2, BLG-CK3, BLG-CK4) 和 接种 根 肿 菌 ‘BLG-S1、BLG-S2、BLG-S3、BLG-S4) 
的 小 叶 获 蓝 根部 ( 根 及 根茎 ) 和 叶子 Can et al., 2019) , 3 个 生物 学 重复 , 并 液 氮 冷冻 , -80 °C 
储存 。 
1.3 方法 
1.3.1 病情 形态 分 级 和 组 织 学 观察 

病情 分 级 参照 农业 部 公益 性 行业 科研 专项 (201003029) 制定 的 标准 ( 杨 华 等 ，2014) 。 
接种 根 肿 菌 后 每 24 h 根据 病情 形态 分 级 观察 根 肿 情况 并 取 小 叶 获 蓝 根 部 进行 组 织 学 观察 ， 
连续 观察 21 d。 将 小 叶 徐 蓝 根部 洗 干 净 ， 取 地 下 1~ 2cm， 于 FAA 固定 液 至 少 固定 24 h; 
后 将 组 织 从 FAA 固定 液 取 出 切 成 薄片 ; 台 盼 蓝 染 液 中 染色 2 min, 染色 后 用 去 离子 水 (ddH20) 
洗 去 残留 的 染 液 ， 清 洗 2~ 3 次 ;， 染 上 色 的 根 组 织 薄 片 放 在 载 玻 片上 进行 显 微 观察 。 
1.3.2 生理 指标 的 测定 

徐 蓝 叶片 生理 指标 的 测定 主要 参照 高 俊 凤 (2006) 和 王 文 龙 (2014) 的 方法 , 具体 如 下 : 
C1) 硫 代 巴 比 妥 酸 显 色 法 检测 MDA; (2) 考 马 斯 亮 蓝 G-250 染色 法 检测 可 溶性 蛋白 ; (3) 
愈 创 木 酚 比 色 法 检测 POD; (4) 氮 蓝 四 唑 光 还 原 法 检测 SOD; (5) 邻 茶 二 酚 法 检测 PPO 
(6) 紫外 分 光 光 度 法 检测 CAT. 
1.3.3 UPLC-MS/MS 测定 
1.3.3.1 样品 制备 

将 小 叶 藉 蓝 的 根 及 根茎 放置 于 冻 干 机 (Scientz-100F〉 中 真空 冷冻 干燥 后 研磨 ， 称 取 50 
mg 的 粉末 , 加 入 70% 甲醇 (1.2 mL), 每 30 min 涡 旋 一 次 (30 s/ 次 , HE 6 次 ); 12 000 rpm, 
离心 3 min 后 ，0.22 pm 微 孔 滤 膜 将 上 清 过 滤 在 进 样 瓶 中 ， 样 品 制备 完成 用 于 后 续 实 验 。 
1.3.3.2 色谱 质谱 条 件 

ExionLC™ AD 超 高 效 液 相 色谱 (https://sciex.com.cn/), Applied Biosystems 4 500 QTRAP 
串联 质谱 Chttps://sciex.com.cn/) 。 液 相 条 件 : 色谱 柱 : SB-C18 (1.8 um, 2.1 mm * 100 mm, 
Agilent); 加 入 了 0.1% 甲 酸 的 超 纯 水 (A) 与 0.1% 甲 酸 乙 且 (B ) 为 流动 相 , 流速 :0.35 mL':min-1; 
洗 脱 梯度 (B 相 比 例 ) : 0.00 min, 5%; 9.00 min 内 线性 增加 至 95%， 并 维持 在 95% 1 min; 
10.00~ 11.10 min, 95%~ 5%; 11.10~ 14 min, 5%; 进 样 量 4 uL; 设 定 柱 温 为 40 °C. MÉ 


in? 


ar 


a3 


条 件 : 电 喷 雾 离 子 源 ， 电 压 : +5500/-4500V; 温度 : 550°C; 离子 源 气体 TI， 气体 OR 
气 帘 气 分 别 设置 50、60 和 25 psi; 碰撞 诱导 电离 参数 : 高。 

OPLS-DA 模型 中 的 VIP 值 结合 差异 倍数 值 fold change 来 筛选 差异 代谢 物 。 筛 选 标准 : 
符合 VIP > 1、fold change > 1.5 和 fold change < 0.67 的 代谢 物 具有 显著 差异 。 


1.3.4 RNA 提取 与 文库 构建 


RNA 提取 与 测序 由 武汉 迈 维 代谢 有 限 公 司 完成 。 采 
RNA Library Prep Kit 试剂 盒 提取 total RNA. 
计 检 测 RNA 完整 性 和 纯度 。 最 终 达标 的 样品 
1.3.5 转录 组 测序 、 分 析 及 注释 

利用 Humina 对 小 叶 获 蓝 根 及 根茎 转录 组 文库 进 
别 将 原始 的 图 像 数 据 转 化 为 原始 数据 Craw data〉， 经 数据 订 
控 获 得 高 质量 序列 (clean reads) , Trinity 组 装 3 


J 高 通 量 测序 。 


| Ilumina 的 NEBNext® UltraTM 
琼脂 糖 凝 胶 电 泳 和 Nano Photometer 分 光 光 度 
| FR EZ) IP SE cDNA 文库 。 


通过 CASAVA 碱 基 识 
人 估 、 过 滤 除 杂 及 宛 余 处 理 等 质 


层次 聚 类 后 得 Unigene。 使 用 DIAMOND 


BLASTX 与 HMMER 软件 将 Unigene 序 列 与 Gene Ontology(GO ), Trembl, Kyoto Encyclopedia 
of Genes and Genomes (KEGG) 、Swiss-Prot、euKaryotic Ortholog Groups (KOG) . NR, 
Protein family (Pfam) 数据 库 进 行 比 对 ， 得 到 基因 功能 信息 。 用 DEGSeq R 包 鉴 定 差异 表达 


基因 


进行 功能 富 集 分 析 。 
1.3.6 qRT-PCR 分 析 


FDR) <0.05, log: 差异 倍数 (fold change, FC) |> 1 的 基因 确定 为 差异 表达 的 基因 3 


(differentially expressed genes, DEGs) ,规定 同时 符合 错误 发 现 率 (false discovery rate, 


对 其 


植物 RNA 提取 同 1.3.4 节 。 根据 转录 组 数据 结果 中 的 FPKM 值 筛选 13 个 差异 基因 ,， 采 


4¥ qRT-PCR. Ail 


分 析 。 扩 增 程序 : 95°C, 2 min 预 变性 ; 
环 ) ; 熔 解 曲线 为 从 65 °C 升 至 95 'C， 每 升温 0.5 C 采集 一 次 荧光 信号 。 每 个 样品 均 设置 
EPCR 结果 的 计算 。 将 实验 结果 数据 表示 为 平均 值 


3 个 重复 , 使 用 2AT 法 进行 实时 荧光 定 引 


用 Monad 试剂 盒 合成 CDNA, LA cDNA 为 模板 、 徐 赣 Actin 基因 (Qu etal., 2019) 为 内 参 进 


Primer-BLAST 设计 引物 ( 表 1)、ABI 7500 荧光 定量 PCR 仪 进行 qRT-PC 
95 'C，5s 变性 ; 60 'C，30s 退火 / 延 人 


(40 个 4 


R 


A 


+ 标准 差 , 并 采用 Graphpad 软件 对 各 组 数据 进行 工 检验 ,以 已 < 0.05 表示 差异 有 统计 学 意义 。 


表 1 
Table 1 


qRT-PCR 引物 
qRT-PCR primers 


HER] ID 基因 名 称 正 向 引物 (5-39 反 向 引物 (5-3) 
Gene ID Gene name Forward primer (5'-3') Reverse primer (5'-3') 
Cluster-16885.9 GH3 TTCGCTTCAGGTATGGTCCG ACCATAGCCACCGAGTTTCG 
Cluster-22704.8  B-ARR GTAAGAGGACTCGGATCGGC CCGTCTCTGCTCTGTTGCAT 
Cluster-22719.0 TGA GAGGGTTTCGTCCATCCGAG GGACTCGTTCACCGCATCAA 
Cluster-27535.0 SAUR ACAACAGCAAACAAGGGATCA GGCAAGGGATCGTGATGGTA 
Cluster-28287.0 AHK TCATCTCCAGCAACGCTCAA CTCTGCACAAACCACTTCCG 
Cluster-37953.0 PP2C CCTGACGTGTCCTATCGACG TTCCATCTCCAACTCCGCAC 
Cluster-35278.2  PP2C CGGTCTTTGGGACGTTGTCT ACGCTTCCTCACACGCTTTA 
Cluster-34981.0 AUX/IAA CTCCGGTTCGTTCGTTCAGA TTCCCTCCCTCTTAGGCTCG 
Cluster-37780.0 | PYR/PYL TCGGTGATCCTGAAATCGGC ACCGAGTATGTGTTCGTCGT 
Cluster-37916.2 ABF GCTTCGGTTCACCAACATCG AGAAGAGAGCCGTGGTGAGT 
Cluster-36994.3 THS CCGAGGCGTCTAGTTTGGTC CAGGAGACGATAGAGCCGAC 
Cluster-36129.3 THS ATATGGCGGTGACGCTCAG CGGTTCGCTTCTCCCATATCA 
Cluster-24362.0 THS ATTGGATGGAAAGCAGGAGGT GTGATCTTGCAGACGCATCC 
内 参 基 基 Actin GCTCACGGAAGCACCT CGACCACTAGCGTAAAGT 


2 结果 与 分 析 
2.1 病情 形态 分 级 与 组 织 学 观察 


根 肿 菌 侵 染 后 植株 的 病情 形态 分 级 和 双 


根部 发 育 正 常 无 肿瘤 ,病情 指数 为 0 级 ， 
接 菌 后 7d 根部 开始 出 ] 
挤 压 变形 (图 1: F, K2: C); 接 菌 后 


昌 织 学 观察 结果 如 图 1 和 图 2 所 示 。 接 菌 后 0 d， 


且 组 织 学 观察 其 细胞 排列 整齐 (图 1: E, 图 2: A); 


岗 肿 大 症状 ， 病 情 指数 为 1 级 ， 组 织 学 观察 表明 皮层 细胞 明显 增 大 ， 


14d 主根 肿块 直径 小 F 葵 基 部 直径 的 2 倍 ， 病情 指 


数 为 3 级 ， 同 时 次 级 游 动 钨 子 遍布 皮层 负 


胞 并 开始 入 侵 维 管 柱 (图 1: G， 图 2: E); 接 菌 
部 直径 的 2~ 3 倍 判定 病情 指数 为 5 级 ， 同 时 次 级 


后 21 d 根部 出 现 较 大 肿瘤 ， 其 直径 是 茎 基 
Ki 


Wah te Tie Ai RE HE (图 1: H, 


2: F) o UeeR RM AO BAW IRAN TARH ea 


侵 染 植株 的 整个 发 病 过 程 ， 同 时 也 说 明 J 
0d 7d 


BLG-CK 


BLG-S 


ARAT E 


人 工 接种 的 可 行 性 。 
14d 


图 1 接 菌 后 0、7、14、21 d 获 蓝 病情 形态 分 级 
Fig. 1 Disease morphological classification of /satis indigotica at 0,7,14 and 21 d after 


inoculation 


A. 接 菌 后 0d RR MAAS, B. 接 菌 后 4d 获 蓝 横 截面 组 织 学 观察 ; C. 接 菌 后 7 d 获 蓝 横 截面 组 
织 学 观察 ; D. 接 菌 后 11 d 获 蓝 横 截面 组 织 学 观察 ; E. 接 菌 后 14 d 获 蓝 横 截 面 组 织 学 观察 ; F. 接 菌 后 21 d 
徐 蓝 横 和 截面 组 织 学 观察 ， 箭头 指向 发 病 部 位 。 


A. Histological observation of cross section of Isatis indigotica on 0 d after inoculation; B. Histological 


observation of cross section of Isatis indigotica on 4 d after inoculation; C. Histological observation of cross 
section of Isatis indigotica on 7 d after inoculation; D. Histological observation of cross section of [satis indigotica 
on 11 d after inoculation; E. Histological observation of cross section of [satis indigotica on 14 d after inoculation; 
F. Histological observation of cross section of 1satis indigotica on 21 d after inoculation. The arrows pointed to the 
site of disease. 
2 WR Daal 152 A Jew PS WS 2 WS 
Fig. 2 Histological observation of Isatis indigotica after infection with Plasmodiophora 


brassicae 


2.2 生理 生化 指标 测定 

由 接 菌 与 未 接 菌 处 理 的 生理 生化 指标 变化 规律 〈 图 3) 可 知 ， 接 菌 处 理 的 可 溶性 和 蛋白、 
MDA 含量 和 防御 酶 SOD. POD, PPO, CAT 活性 随 接 菌 时 间 的 延长 呈 不 断 上 升 的 趋势 ; 
CAT 和 SOD 活性 在 接 菌 后 7~ 21 d， 实 验 组 均 显著 高 于 相应 对 照 组 (P< 0.05)〉 ; 可 溶性 重 
白 \ 丙 二 醛 含量 以 及 POD、PPO 活性 在 接 菌 后 14~ 21 d, 实验 组 均 显 著 高 于 对 照 组 (P < 0.05). 
以 上 结果 说 明 获 蓝 为 抵抗 根 肿 菌 侵 染 自身 做 出 了 防御 反应 。 
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*P < 0.05, “P< 0.01, P< 0.001 manifest that there is a noteworthy discrepancy between the experimental group 


and the control group at the same time point. 


3 ”接种 根 肿 黄 


后 


不 同时 间 松 监 叶 生 理 生 化 指标 测定 


Fig. 3 Determination of physiological and biochemical indexes of Isatis indigotica at different 


time after inoculation with Plasmodiophora brassicae 


2.3 代谢 组 分 析 
2.3.1 样品 检测 及 OPLS-DA 分 析 


本 实验 对 7 组 样本 CBLG-CK1 (S1) 、BLG-CK2、BLG-CK3、BLG-CK4、BLG-S2、 


BLG-S3、BLG-S4) 进行 代谢 研究 。 基 于 检测 习 


EF 台 和 自 建 数 据 库 共 检 测 到 161 种 生物 碱 ， 其 


他 类 生物 碱 数量 最 多 ， 占 总 生物 碱 的 53%， 其 次 是 嘻 吹 类 生物 碱 占 总 生物 碱 的 32% (图 4)。 


OPLS-DA 正 交 偏 最 小 二 乘 判 别 分 析 (Orthogonal partial least squares-discriminant analysis ) 结 
果 ( 图 5) 显示 各 组 样品 均 能 得 到 明显 分 离 ， 且 BLG-CK3 和 BLG-S3 的 分 离 吕 
BLG-CK4 和 BLG-S4 次 之 ， 这 说 明 接 菌 


53% 


14 d 时 根 肿 菌 对 获 蓝 代谢 影响 最 大 。 
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图 4 代谢 物 类 别 分 析 


Fig.4 Class analysis of metabolites 


E 离 最 大 ， 
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图 5 接种 根 肿 菌 后 7、14、21 d 及 其 相应 对 照 组 的 OPLS-DA 得 分 图 
Fig.3 OPLS-DA score map at 7,14, 21 d after inoculation with Plasmodiophora brassicae and 


corresponding control group 


2.3.2 差异 代谢 物 筛 ; 

进一步 对 接 菌 后 7、14、21 d 的 获 蓝 差异 代谢 物 进行 分 析 ， 结 果 发 现 : 
BLG-CK2-vs-BLG-S2 有 16 种 差异 代谢 物 ，6 个 上 调 ，10 个 下 调 ，BLG-CK3-vs-BLG-S3 
17 种 差异 代谢 物 ，9 个 上 调 ，8 个 下 调 ; BLG-CK4-vs-BLG-S4 有 39 种 差异 代谢 物 ，32 个 上 
调 ,7 个 下 调 ;BLG-S2-vs-BLG-S3 有 19 个 差异 代谢 物 ,3 个 上 调 ,16 个 下 调 ; BG-S2-vs-BLG-S4 
有 45 个 差异 代谢 物 ，19 个 上 调 ，26 个 下 调 ; BG-S3-vs-BLG-S4 有 34 个 差异 代谢 物 ，14 个 
上 调 ，20 个 下 调 〈 见 表 2) 。 如 图 6 所 示 , 将 BLG-CK2-vs-BLG-S2、BLG-CK3-vs-BLG-S3、 
BLG-CK4-vs-BLG-S4 这 3 组 差异 代谢 物 取 交集 发 现 共 有 5 个 差异 代谢 物 , DWNT 
Isatindosulfonic acid B、5，6- 二 羟基 叮 唆 -$-O-B- 葡 萄 糖苷 、 泛 酰 统 基 乙 胶 和 对 香 豆 酰 亚 精 腕 ， 
其 中 前 3 种 为 呀 唆 类 生物 碱 . 将 BLG-S2-vs-BLG-S3、BLG-S2-vs-BLG-S4、BLG-S3-vs-BLG-S4 
这 3 组 差异 代谢 物 取 交集 发 现 共 有 2 个 差异 代谢 物 , 分 别 是 呀 唆 类 生物 碱 Isatisindigoticanine 
B 和 唆 啉 类 生物 碱 2- 氧 -3，4- 二 氧 -1 H- 唆 啉 -3- 羧 酸 。 


表 2 差异 代谢 物 数量 统计 


Table 2 Differential metabolite count statistics 


比 对 组 总 数 上 调 下 调 

Group Total Up Down 
BLG-CK2-vs-BLG-S2 16 6 10 
BLG-CK3-vs-BLG-S3 17 9 8 
BLG-CK4-vs-BLG-S4 39 32 T 
BLG-S2-vs-BLG-S3 19 3 16 
BLG-S2-vs-BLG-S4 45 19 26 


BLG-S3-vs-BLG-S4 34 14 20 


BLG-S2-vs-BLG-S3 BLG-S2-vs-BLG-S4 


BLG-CK3-vs-BLG-S3 


BLG-CK2-vs-BLG-S2 


BLG-CK4-vs-BLG-S4 BLG-S3-vs-BLG-S4 


6 差异 代谢 物 韦 恩 图 
Fig.6 Venn diagram of differential metabolite 


2.3.3 差异 代谢 物 KEGG 富 集 分 析 

KEGG 通路 富 集 分 析 结 果 见 图 7。 结 果 表 明 , 各 组 别 富 集 在 代谢 通路 (Metabolic pathways, 
ko001100) 上 的 差异 代谢 物 最 多 ， 且 生物 碱 和 氨基酸 相关 的 代谢 通路 变化 明显 ; 
BLG-CK2-vs-BLG-S2、BLG-CK3-vs-BLG-S3、BLG-CK4-vs-BLG-S4 这 3 组 都 富 集 在 托 烷 、 
哌 啶 和 吡啶 生物 碱 生物 合 成 途径 (Tropane，piperidine and pyridine alkaloid biosynthesis , 
ko00960) ; 除了 BLG-CK3-vs-BLG-S3, 其 他 组 都 在 色 和 氨 酸 代谢 途径 (Tryptophan metabolism, 


ko003800) 上 富 集 ， 色 氨 酸 合成 途径 是 叫 吹 类 生物 碱 合成 的 前 体 途 径 。 
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差异 代谢 物 KEGG 代谢 通路 富 集散 点 分 析 图 


Fig.7 Scatter map of differential metabolites KEGG pathway enrichment 


2.4 RNA-Seq 分 析 
2.4.1 转录 组 数据 组 装 与 质量 分 析 

为 了 进一步 挖掘 根 肿 菌 胁迫 下 黎 蓝 生物 碱 积累 机 制 ,对 徐 蓝 21 个 样本 进行 转录 组 测序 
和 分 析 得 到 1 200 974 786 条 raw reads 和 1 171 808 210 个 clean reads， 共 获得 175.77 Gb 的 
效 数 据 ;， 各 样本 Clean Data 均 达 到 7 Gb, Q20 (Qphred 值 不 低 于 20 的 碱 基数 占 总 碱 基数 
的 百分比 ) 碱 基 百 分 比 均 在 97% 以 上 ，Q30 (Qphred 值 不 低 于 30 的 碱 基 数 占 总 碱 基数 的 百 
分 比 ) 碱 基 百 分 比 均 在 92% 以 上 ;GC 舍 量 (G 和 C 占 总 碱 基 数量 的 百分比 ) 是 47.0%~ 47.69% 
Trinity 拼接 得 到 83 775 个 Unigene， 平 均 长 度 为 1 708 bp， 最 长 达到 16 523 bp， 最 短 是 201 
bp, N50 是 2367 bp. Unigene 长 度 分 布 图 (图 8) 显示 ，51 799 条 (61.83%) Unigene KE 
超过 1000 bp, 25 874 条 (30.89%) Unigene 长 度 超 过 2 000 bp。 以 上 说 明 转录 组 数据 质量 
较 高 ， 可 进行 后 续 分 析 。 
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8 徐 蓝 转录 组 Unigene 长 度 分 布 图 
Fig. 8 Length distribution of Unigene in transcriptome of Isatis indigotica 


数量 Number 
S 
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2.4.2 Unigene 的 功能 注释 
用 NR、TrEMBL、GO、Swiss-port、KOG、KEGG、Pfam 共 7 个 数据 库 进 行 Unigene 
注释 , 分 别 成 功 注释 69 350 (82.78% ) ~ 70 281 (83.89%) ~ 60 538 (72.26%) ~ 54 696 (65.29%). 
44 620 (53.26%) . 55271 (65.98%) . 55387 (66.11%) 个 基因 ， 至 少 在 一 个 数据 库 注释 
成 功 的 Unigene 有 72 965 条 (87.1%) 。 而 样本 间 基 因 表达 水 平 的 相关 性 分 析 表 明 样 本 在 一 
组 中 具有 较 高 的 一 致 性 ， 确保 了 后 面 分 析 结 吉 果 的 可 靠 性 。 
在 本 研究 中 ， 共 有 60538 条 Unigene 得 到 了 GO 注释 。 总 体 分 为 3 KA: 分 子 功能 
(molecular function ) 中 主要 富 集 条 目 有 催化 活性 (catalytic activity, 30 758 A) AZ FA (binding, 
36 084 条 ) ; 细胞 组 分 (cellular component〉 中 细胞 解剖 实体 (cellular anatomical entity ) 
Unigene 数量 最 多 ，51 807 条 ; 生物 过 程 (biological process) 中 Unigene 富 集 较 多 的 类 别 有 : 
细胞 过 程 (cellular process，40 504 条 ) 、 应 激 适 应 (response to stimulus, 18 474) 、 代 谢 
过 程 (metabolic process, 32 229 条 ) 、 生 物 调控 (biological regulation, 16 606 条 ) (图 9)。 
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1. 细胞 过 程 (40504) ; 2. 代谢 过 程 (32229) ; 3. 刺激 响应 (18474) ; 4. 生物 调节 (16 606) ; 5. Æ 
物 过 程 的 调节 (14988) ; 6. 发 展 过 程 (9177) ; 7. 多 细胞 生物 过 程 (8 714) ; 8. 定位 (7499) ; 9. 信 
号 (5694) ; 10. 复制 (5 445) ; 11. 生殖 过 程 (5415) ; 12. 种 间 相 互 作 用 的 生物 过 程 (4 780) ; 13. Æ 
物 过 程 的 负 调控 (3 834) ; 14. 积极 调节 生物 过 程 (3 713) ; 15. 生长 (2 032) ; 16. 免疫 系统 过 程 (1 382); 
17. 多 生物 过 程 (1 319) ; 18. 节律 过 程 (570) ; 19. 解毒 (232) ; 20. 运动 (162) ; 21. 生物 黏附 〈100); 
22. 种 内 相互 作用 的 生物 过 程 (64) ; 23. 氮 利 用 率 (46) ; 24. WAHR (35) ; 25. 色素 沉积 〈14) ; 

26. 硫 利用 率 〈1) ; 27. 碳水 化 合 物 利 用 率 (1) ; 28. 细胞 解剖 实体 (51 807) ; 29. 包含 蛋白 复合 物 (9 
458) ; 30. 结合 (36 084) ; 31. 催化 活性 (30 758) ; 32. 转运 体 活性 (4 183) ; 33. 转录 调节 活性 (3 266); 
34. 结构 分 子 活性 (2 128) ; 35. 分 子 功能 调节 (1747) ; 36. 分 子 转 导 活性 (748) ; 37. 翻译 调节 活性 
(639) ; 38. 抗 氧化 活性 (470)〉; 39. 分 子 载体 活性 (213〉 ; 40. 分 子 衔接 活性 (209) ; 41. RAS 
FHB (121) ; 42. 和 蛋白质 标签 (76) ; 43. 一 般 起 始 因子 活性 〈59) ; 44. 小 分 子 传感器 活性 (42) ; 45. 
养 库 活动 (40) ; 46. 毒素 活性 (5) 。 

1. Cellular process (40 504); 2. Metabolic process (32 229); 3. Response to stimulus (18 474); 4. Biological 


[ny 


m 


regulation (16 606); 5. Regulation of biological process (14 988); 6. Developmental process (9 177); 7. 
Multicellular organismal process (8 714); 8. Localization (7 499); 9. Signaling (5 694); 10. Reproduction (5 445); 
11. Reproductive process (5 415); 12. Biological process involved in interspecies interaction between organisms (4 
780); 13. Negative regulation of biological process (3 834); 14. Positive regulation of biological process (3 713); 
15. Growth (2 032); 16. Immune system process (1 382); 17. Multi-organism process (1 319); 18. Rhythmic 
process (570); 19. Detoxification (232); 20. Locomotion (162); 21. Biological adhesion (100); 22. Biological 
process involved in intraspecies interaction between organisms (64); 23. Nitrogen utilization (46); 24. Carbon 
utilization (35); 25. Pigmentation (14); 26. Sulfur utilization (1); 27. Carbohydrate utilization (1); 28. Cellular 
anatomical entity (51 807); 29. Protein-containing complex (9 458); 30. Binding (36 084); 31. Catalytic activity 
(30 758); 32. Transporter activity (4 183); 33. Transcription regulator activity (3 266); 34. Structural molecule 
activity (2 128); 35. Molecular function regulator (1 747); 36. Molecular transducer activity (748); 37. Translation 
regulator activity (639); 38. Antioxidant activity (470); 39. Molecular carrier activity (213); 40. Molecular adaptor 
activity (209); 41. Protein folding chaperone (121); 42. Protein tag (76); 43. General transcription initiation factor 
activity (59); 44. Small molecule sensor activity (42); 45. Nutrient reservoir activity (40); 46. Toxin activity (5). 
9 Unigene 的 GO 注释 
Fig.9 GO annotation distribution of Unigene 


KOG 功能 分 类 分 析 《〈 图 10) 可 进一步 了 解 获 蓝 转 录 组 的 功能 ， 本 次 共有 25 个 不 同 功 
能 且 种 类 齐全 类 群 , 包含 大 多 数 的 生命 活动 。 其 中 10 102 条 被 富 集 到 一 般 功 能 预测 , 是 KOG 


类 群 中 Unigene 数量 最 多 的 ; FPA ONE a Eth. AeA SHA, 4932 条 ; 
注释 到 信号 转 导 功能 的 也 有 4 461 个 基因 ， 其 他 功能 类 群 的 基因 丰 度 也 不 尽 相 同 。 


10 000 


7 500 


基因 数量 Number of genes 


0 
ABCDEFGHIJKLMNOPQRSTUVWYZ 
分 类 Classification 


A. RNA 加 工 和 修饰 (2 023) ; B. 染色 体 结构 和 活力 (715) ; C. 能 量 生 成 和 转换 (1 875) ; D. 细胞 周 
期 控制 、 蛋 白质 折 县 和 分 子 伴侣 (1 269) ; E 氨基 酸 运输 和 代谢 (1 905); F. 核 苷 酸 运输 和 代谢 (475); 
G. 碳水 化 合 物 运输 和 代谢 (2335); H. 辅酶 运输 和 代谢 〈520) ; L 脂 质 运输 和 代谢 《1 797) ; J. 翻译 、 
核糖 体 结构 和 生物 合成 (3 214) ; K. 转录 (2686); L. 重复 、 重 组 和 修饰 (1 643); M. 细胞 壁 膜 生 物 
合成 642) ; N. 细胞 运动 (24) ; O. 翻译 后 修饰 、 蛋 白质 折 肝 和 分 子 伴侣 (4 932)〉 ; P. 矿 脂 运输 和 代 
谢 (1288) ; Q. 次 生 代 谢 物 合成 、 运 输 和 代谢 (1 558); R. 一 般 功能 预测 (10 102) ; S. 功能 未 知 (2 450); 
T. 信号 传导 (4461); U. 细胞 内 转运 、 分 泌 和 小 泡 运 输 (2096) ; V. 防卫 机 制 (357) ; W. 胞 外 结构 
(104) ; Y. 核 结 构 (138) ; Z. 细胞 构架 (1135) 。 

A. RNA processing and modification (2 023); B. Chromatin structure and dynamics (715); C. Energy production 


EN 


and conversion (1 875); D. Cell cycle control, cell division, chromosome partitioning (1 269); E. Amino acid 
transport and metabolism (1 905); F. Nucleotide transport and metabolism (475); G. Carbohydrate transport and 
metabolism (2 335); H. Coenzyme transport and metabolism (520); I. Lipid transport and metabolism (1 897); J. 
Translation, ribosomal structure and biogenesis (3 214); K. Transcription (2 686); L. Replication, recombination 
and repair (1 643); M. Cell wall/membrane/envelope biogenesis (642); N. Cell motility (24); O. Posttranslational 
modification, protein turnover, chaperones (4 932); P. Inorganic ion transport and metabolism (1 288); Q. 
Secondary metabolites biosynthesis, transport and catabolism (1 559); R. General function prediction only (10 
102); S. Function unknown (2 450); T. Signal transduction mechanisms (4 461); U. Intracellular trafficking, 
secretion, and vesicular transport (2 096); V. Defense mechanisms (357); W. Extracellular structures (104); Y. 
Nuclear structure (138); Z. Cytoskeleton (1 135). 

图 10 Unigene 的 KOG 分 类 

Fig. 10 KOG classification of Unigene 


2.4.3 差异 表达 基因 的 筛选 与 功能 分 析 

PS eA 7. 14. 21 d 与 同时 间 对 照 组 比较 分 别 有 2 439、256、6 437 个 差异 表达 基 
因 ， 结 果 表 明和 其 蓝 在 根 肿 菌 侵 染 后 基因 表达 发 生 显著 变化 ， 其 中 21 d 的 DEGs 最 多 ， 说 明 
根 肿 菌 侵 染 后 期 获 蓝 基因 表达 发 生 剧 烈 变化 ; 这 3 组 共有 17 个 差异 表达 基因 ， 其 中 9 个 是 
上 调 的 ，8 个 下 调 。 接 菌 后 7、14、21 d 的 实验 组 随 着 接 菌 时 间 的 延长 ， 差 异 基因 数目 在 不 
断 增 加 ， 其 对 根 肿 菌 的 响应 更 加 剧烈 〈 表 3) 。 


<< 


表 3 差异 基因 数量 统计 


Table 3 Differentially expressed genes in different samples 


比 对 组 总 数 上 调 下 调 

Group Total Up Down 
BLG-CK2-vs-BLG-S2 2 439 1 298 1141 
BLG-CK3-vs-BLG-S3 256 116 95 
BLG-CK4-vs-BLG-S4 6 437 5 838 599 
BLG-S2-vs-BLG-S3 2 451 1174 1277 
BLG-S2-vs-BLG-S4 8 752 6 494 2 256 
BLG-S3-vs-BLG-S4 6 213 5 340 853 


对 DEGs 进行 GO term 富 集 分 析 , 根据 细胞 组 成 分 子 功能 和 生物 代谢 对 基因 进行 分 类 ， 
EA 9 805 个 DEGs 被 富 集 到 这 3 个 GO 类 别 中 。 细 胞 组 成 过 程 中 细胞 解剖 实体 〈cellular 
anatomical entity) 的 DEGs 居 第 一 位 置 ， 细胞 过 程 〈cellular process) 、 刺 激 反 应 (response 
to stimulus) 和 代谢 过 程 (metabolic process) 是 DEGs 在 生物 过 程 中 富 集 较 多 的 3 个 功能 ; 
结合 (binding) 和 催化 活性 Catalytic activity) 是 DEGs 在 分 子 功能 中 的 主要 功能 ( 表 4) 。 


表 4 差异 基因 数 及 功能 注释 


Table 4 Function annotation and gene dosage of the differentially expression genes 


分 类 功能 注释 总 数 比 对 组 


Classification Function annotation Total Contrast group 


BLG-CK2- BLG-CK3- BLG-CK4- BLG-S2- BLG-S2- BLG-S3- 
VS- VS- VS- VS- VS- VS- 


BLG-S2 BLG-S3 BLG-S4 BLG-S3 BLG-S4 BLG-S4 


细胞 组 成 细胞 解剖 实体 16 973 1871 181 3 821 1 845 5 584 3 671 
Cellular Cellular anatomical 


component entity 


包含 蛋白 复合 物 3 488 168 18 1 042 131 1 137 992 


Protein-containing 
complex 
生物 过 程 细胞 过 程 13 688 1 415 140 3 244 1 338 4 463 3 088 
Biological Cellular process 
process 代谢 过 程 10 692 1 107 105 2 561 1 034 3 482 2 403 
Metabolic process 
刺激 反应 6 495 852 93 1 301 850 2 139 1 260 
Response to 


stimulus 


生物 调节 5 621 693 63 1 169 694 1 856 1 146 
Biological 


regulation 


并 


生物 过 程 调控 5 054 629 57 1 046 628 1 672 1 022 
Regulation of 


biological process 


分 子 功能 
Molecular 


function 


由 于 徐 蓝 中 的 主要 活性 成 分 为 生物 碱 , 氨基 酸 又 是 生物 碱 合成 的 基 而 
酸 和 生物 碱 相关 的 代谢 通路 。 为 了 更 清 
E, X} BLG-CK2-vs-BLG-S2、BLG-CK3-vs-BLG-S3、BLG-CK4-vs-BLG-S4 
进行 KEGG 富 集 分 析 ， 结 果 见 表 5。 结 果 发 现 3 组 共同 富 集 在 11 个 生物 碱 
其 中 色 氨 酸 代谢 途径 富 集 的 差异 基 


通路 上 存在 相关 改 
这 3 组 差异 基因 
相关 的 代谢 通路 上 ， 


结合 

Binding 
催化 活性 
Atalytic activity 


转运 蛋白 活性 


Transporter activity 


翻译 调节 蛋白 活性 


Transcription 


regulator activity 


11 912 
10 329 
1 348 163 
1 058 187 


1300 


1113 


121 2 703 1 237 
166 2315 1121 
23 265 172 
14 160 173 


因数 量 较 多 。 


3931 2 620 
3395 2219 
463 262 
370 154 


Hy 故 主要 分 析 氨 基 


楚 地 了 解 所 关注 的 生物 碱 类 化 合 物 与 差异 基因 是 否 在 


T 


表 5 生物 碱 合成 相关 通路 KEGG ER 
Table 5 KEGG enrichment in pathways related to alkaloid synthesis 
Ko-ID KEGG 通路 名 称 差异 表达 基因 数量 
KEGG pathway name DEGs count 
BLG-CK2-vs- | BLG-CK3-vs- | BLG-CK4-vs- 
BLG-S2 BLG-S3 BLG-S4 

ko00950 — FRMEENRKAE Py tate YA 5 2 8 
Isoquinoline alkaloid biosynthesis 

ko00960 Eki, PRUEIME AE Ee BM 4 1 8 
Tropane ， piperidine and pyridine alkaloid 
biosynthesis 

ko00360 ARAN ART 5 1 8 
Phenylalanine metabolism 

ko00350 ” 栈 氨 酸 代谢 6 2 18 
Tyrosine metabolism 

ko00380 GARRI 6 2 25 
Tryptophan metabolism 

ko00300 ” 赖 氨 酸 生 物 合成 1 0 5 
Lysine biosynthesis 

ko00330 ” 精 氨 酸 和 且 氨 酸 代谢 14 1 24 
Arginine and proline metabolism 

ko00400 NAR MAR EAREWAR 10 1 18 
Phenylalanine, tyrosine and tryptophan 
biosynthesis 

ko00220 ” 精 氨 酸 生 物 合成 11 1 18 
Arginine biosynthesis 

ko00901 ” 呵 哄 生物 碱 生物 合成 0 0 3 
Indole alkaloid biosynthesis 

ko00996 ”各 种 生物 碱 生物 合成 0 0 1 


Biosynthesis of various alkaloids 


2.5 生物 碱 合成 途径 相关 基因 挖掘 


本 研究 重点 对 呀 唆 类 生物 碱 合 成 途径 及 相关 基因 进行 挖掘 。 结 合 图 11 分 析 : 分 支 酸 可 
在 邻 氨基 葵 甲 酸 合 酶 (Anthranilate synthase, AS)、 色 和 氨 酸 合成 酶 (Trytophan synthase asubunit, 
TSA) 等 酶 的 反应 下 合成 叫 唆 ; 最 后 经 过 一 系列 反应 合成 色 氨 酸 。 色 氨 酸 脱羧 酶 (L-tryptophan 


decarboxylase, TDC) 是 至 关 重 要 的 酶 ， 色 氢 酸 可 在 其 催化 下 生成 色 胶 ， 并 在 叫 唆 -3- 丙 酮 酸 
单 加 氧 酶 (indole-3-pyruvate monooxygenase, YUCCA) 的 催化 作用 下 生成 呀 唆 -3- 乙 酸 酯 ; 
fe tH Hy HE PRE Hh, SUE (aldehyde dehydrogenase, ALDH) 的 催化 下 合成 $- 羟 基 呵 吹 -3- 乙 酸 。 


分 支 酸 Chorismate 


AS 


邻 氨 基 茶 甲酸 Anthranilate 


| TSA 


y 
Hg|BE Indole 


L-4 ARA L-tryptophan 


TDC 


_ YUCCA D 
色 胺 Tryptamine 一 一 一 > 日 吹 -3- 乙 酸 酯 Indole-3-acetate 
1 
ICYP71P!1 
5-FRFEN| WE CHE 5-hydroxyindoleacetaldehyde 


ALDH 


S~FRAEMG| Wk CERNE 
ASMT 


5-hydroxyindoleacetate 


5-FA ENS WE LARAS 5-methoxyindoleacetate 


图 11 
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Fig. 11 Biosynthetic pathway of indole alkaloids 


本 研究 从 所 有 的 DEGs 中 筛选 到 18 个 基因 ， 分 别 编码 呵 哄 生物 碱 合 成 途径 的 5 种 关键 
酶 (AS. TSA, TDC, YUCCA, ALDH) (图 12: A) ， 另 外 18 个 基因 编码 异 唆 啉 类 生物 


碱 合 成 途径 的 2 种 关键 酶 ， 带 巴 因 
(tyrosine aminotransferase, TAT) 
表达 量 均 上 调 的 一 共有 9 个 DEGs: 


合成 酶 (thebaine synthase, THS) 和 栈 氨 酸 氮 基 转 移 酶 
(图 12: B) 。 与 未 接 菌 相 比 ， 接 菌 后 7、14、21 d 基因 
3 个 编码 THS 的 基因 Cluster-24362.0. Cluster-36994.3 ~ 


Cluster-36129.3; 2 个 编码 TAT 的 基因 Cluster-31730.0. Cluster-28040.6; 2 个 编码 YUCCA 
的 基因 Cluster-36192.0, Cluster-36192.1; 以 及 2 个 编码 ALDH 的 基因 Cluster-32381.0. 


Cluster-28395.0. 


Cluster-24362.0 K : 
Cluster-36994.3 
Cluster-36129.3 E ° 
Cluster-36269.0 
Cluster-36994.0 i 
Cluster-24922.0 ÎTI 
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Cjusrer-38070.4 2 


Cluster-39788.9 
Cluster-9074.0 行 分 组 


Row_Group Row_Group 
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Cluster-36192.1 ppc Cluster-243 62.3 
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Fj Cluster-32381.0 E LDH Cluster-28040.6 
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Cluster-8907.0 

Cluster-10828.0 
Cluster-13450.0 
Cluster-39075.0 


Cluster-39394.4 
Cluster-10426.0 
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ta im Ke 
Ra as an o; Row_Group 


ERRE A HE m P ARKAE FPKM 标准 化 处 理 后 的 值 。 
The color represents the gene expression in the sample after FPKM normalization. 
12 ”生物 碱 合 成 相关 基因 表达 量 分 析 


Fig. 12 Analysis of gene expression associated to alkaloid synthesis 


2.6 差异 表达 基因 的 qRT-PCR 验证 
为 了 验证 转录 组 数据 的 准确 性 ， 选 择 了 THS (Cluster-24362.0、Cluster-36994.3、 
Cluster-36129.3 ) , AUX/IAA (Cluster-34981.0) , GH3 ( Cluster-16885.9) , SAUR 
(Cluster-27535.0), PYR/PYL (Cluster-37780.0), PP2C (Cluster-35278.2、 Cluster-37953.0), 
ABF (Cluster-37916.2) , B-ARR (Cluster-22704.8) , TGA (Cluster-22719.0) , AHK 
C Cluster-28287.0) 进行 qRT-PCR 验证 (图 13: A) 。 结 果 表 明 ， 上 述 酶 基因 在 
BLG-CK2-vs-BLG-S2、BLG-CK3-vs-BLG-S3、BLG-CK4-vs-BLG-S4 中 的 表达 趋势 与 RNA-seq 
结果 一 致 (P<0.05) (图 13: B) ， 证 明 测序 结果 真实 可 靠 。 
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A. 采用 25T 法 计算 13 个 差异 表达 基因 的 相对 表达 量 ; 


mu 


qRT-PCR results and RNA-seq data. 
图 13 差异 表达 基因 
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B. qRT-PCR 结果 与 RNA-seq 数据 之 间 的 相关 性 。 


A. The relative expression levels of 13 DEGs were calculated by the 244°T method; B. Correlation between 


的 qRT-PCR 验证 
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Fig. 13 qRT-PCR validation of differentially expressed genes 


3 讨论 与 结论 


3.1 根 肿 菌 侵 染 对 黎 蓝 根部 的 影响 


I 


大 量 研究 表明 ， 感 染 根 肿 菌 的 植株 其 根部 会 不 断 膨大 ， 直 至 变 成 大 小 形状 不 一 的 肿瘤 ， 


并 且 根 部 细胞 会 因为 次 级 游 动 钨 子 的 入 
植株 根部 的 肿 大 程度 能 反映 其 病情 指数 , 根部 细 


本 研究 连续 21 d 对 根 肿 菌 侵 染 获 蓝 的 过 程 进行 实时 动态 观察 ,而 
站 数 不 断 增加 , 分 别 为 0 级 、1 级 、3 级 、5 级 ; 组织 学 观察 发 现在 接种 后 7~ 14 


21 d MATS 


d {RI R REPRE RA HL» SR UE, RSE (2022) ANE 


ZT RABI. ARPA (Wei etal., 2021) 。 通 常 


胞 的 变化 情况 能 反映 根 肿 菌 的 整个 侵 染 过 程 。 


究 发 现 接 菌 后 0、7、14、 


完结 果 一 致 ; 


虽然 Wei (2021) 的 研究 在 根 肿 菌 


所 


不 仅 有 利于 从 整体 全 面 的 观察 根 肿 


且 为 后 期 生化 指标 、 转 录 组 
3.2 黎 蓝 响应 根 肿 菌 侵 染 生理 生化 指标 的 变化 


ARH ead 152 


RP, 


PS WEEP 
(Wang et al., 2020) . A 


学 和 代谢 组 


结果 发 现 ， 可 溶性 重 白 和 丙 二 醛 含量 在 根 肿 菌 
侵 染 后 期 ,植物 受到 的 
POD. CAT. PPO 活性 
肿 菌 对 植物 的 伤 
性 于 14 d 出 现 显著 的 差异 , 说 明了 POD 和 PPO 发 挥 作 月 


加 ， 说 明 根 肿 


ics] 


ae 


Fo» 


身 会 做 出 


Pe rahe 


ian 


Nb Fi 
程度 更 严 


究 的 寄主 与 侵 染 条 件 不 同 有 关 。 此 外 , 本 研究 采取 病情 形 


ob H 


侵 染 时 间 上 与 本 研究 结果 有 所 出 入 ， 考 虑 到 这 可 
态 分 级 结合 组 织 学 观察 的 方法 ， 


HG AE 


a URAL, ARATE 
学 取样 时 间 提 供 了 到 


他 专 性 寄生 病原 菌 
论 依据 。 


做 参照 , 而 


mh 


E 


系列 的 应 激 反应 ,这 与 大 多 数学 者 看 
究 通过 对 接 菌 后 0、7、14、21 d 获 蓝 的 生理 生化 指标 进行 测定 ， 
Ua 14d 和 21 d 与 未 接 菌 相 比 
H 


; 根 


完结 果 一 臻 


含量 显著 增 


处 理 下 


松 赣 抗 氧 化 酶 SOD. 


变化 趋势 都 是 不 断 上 升 ， 说 明 徐 蓝 通过 增加 抗 氧化 酶 活性 来 抵御 根 


CAT 和 SOD 活性 在 接 菌 后 7 d 开始 存在 明显 差异 ， 但 POD 和 PPO 活 


来 说 存在 灌 后 性 ， 这 与 郑 光 明 等 (1999) 的 相关 报道 一 致 。 


3.3 根 肿 菌 侵 染 对 获 蓝 生物 碱 合 成 的 影响 


徐 蓝 作为 我 国 的 大 宗 药材 ， 其 吓 唆 类 生物 胡 


E 


成 表 告 依 


HAR ELAR TAI 2 种 防御 酶 SOD. CAT 


aN 


有 丰富 的 药理 活性 。 


长 春花 


(Catharanthus roseus) 〈 林 颖 等 ，2020) ~ BEA (Rauvolfi serpentina) (Dey et al., 2022). 


H (Uncaria rhynchophylla) 


模式 植物 ， 徐 蓝 中 吧 吹 生物 碱 代谢 途 


体 途 径 (Huang et al., 2016) 。 H 
因 功 能 的 研究 少见 报道 , 并且 有 看 
2016; Jazayeri et al., 2022) 。 


分 子 机 制 。 


J, 


( 刘 扬 等 ， 
径 复 杂 多 样 ， 


2021) 


SAY LE AIG RR AE WMA BN 


其 中 


RY 


代谢 组 学 


通 


发 现 不 管 是 接 


其 是 色 和 氨 酸 代谢 , 说 明 根 肿 菌 侵 染 影 
BLG-CK3-vs-BLG-S3、BLG-CK4-vs-BLG-S4 这 3 个 组 别 的 差异 表达 基 
第 选 出 了 11 条 与 生物 碱 合 成 相关 的 代谢 通路 ， 
上 转录 加 代谢 联合 分 析 共 挖掘 18 个 差异 基因 参与 叫 唆 生物 碱 合 成 途 
TDC, YUCCA, ALDH, # 
Lig ital, AE ia. W 
(Cluster-36192.0. Cluster-36192.1) 表达 量 均 升 
E 物 胁迫 的 耐 受 性 ， 但 随 着 根 肿 菌 


路 富 集 分 析 ， 
因数 量 较 多 。 


径 上 游 的 5 种 关键 酶 AS. TSA, 
键 酶 的 位 置 上 存在 差异 表达 基因 


YUCCA 基因 
完结 果 一 致 ， 


ics] 


在 此 基 耐 


过 对 接 菌 
处 理 ( 接 


都 呈现 出 对 


a 


因 的 表达 量 


完 表 明 逆 境 胁 迫 


Sb BZ 
HE RZ 


HIRE HY 


EAN 


骏 代 谢 是 叫 吹 类 生物 碱 合 成 的 前 
成 的 合成 机 制 、 关 键 酶 、 3 
响 活性 成 分 生物 碱 的 积累 ( 


完 的 


多 种 基 
Ate 


青 晓 清 于 > 


茵 和 不 接 


旺 先 升 高 再 降低 ， 这 可 
( Cluster-32381.0, Cluster-28395.0) 基 


Al) 还 是 ] 


因此 ， 有 必要 深入 研究 根 肿 


茵 胁迫 


引 吹 类 生物 碱 合成 的 


PAE 


前 后 7、14、21 d 获 蓝 根 次 级 代谢 产物 生物 碱 进 行 OPLS-DA 分 析 ， 
取样 时 间 不 同 ， 代 谢 物 都 存在 差异 ， 且 次 级 代谢 
物 以 呀 唆 类 生物 碱 居多 。 差 异 代谢 物 KEGG 富 集 分 析 表 明 各 组 多 富 集 在 氨基 酸 代谢 途径 万 
响 生物 碱 类 化 合 物 合 成 。 进 一 步 对 BLG-CK2-vs-BLG-S2、 


因 进 行 KEGG 代谢 通 


中 色 氨 酸 合成 相关 的 差异 基 


YUCCA 和 ALDH 这 2 个 关 


tA ACH GAR BT 


相 比 ， 接 菌 后 


能 是 植物 


因 在 接 


为 醛 脱 氧 酶 基因 


路 。 


Fab, BANE RH TY 


中 被 注释 。 


Cluster-36994.3. Cluster-36129.3 ) 在 根 肿 菌 
证 中 也 具有 相同 的 表达 趋 


期 可 通过 全 基 


i 


da 


aN 2 PES BEM 
目前 THS 酶 的 相关 报 
一 (Ozber et al., 2023) 。 本 研究 结果 显 


较 少 


势 ， 说 明 THS 


高 ， 这 与 Cao (2019) 的 研 
受 染 时 间 的 不 断 增加 ，YUCC4 
身 对 道 境 胁 迫 做 出 


的 应 激 反 应 ; ALDH 


A Ja AHX EEJ 
家 族 可 以 氧化 有 毒 醛 类 物质 ， 降 低 脂 质 过 氧化 ， 提 高 植 4 
(Du et al., 2022; Zhang et al., 2023) 。 由 于 缺乏 根 肿 菌 
唆 类 生物 碱 合 成 通路 中 的 全 部 基因 
究 造 成 困难 。 为 此, 后 


进行 注释 分 析 , BOTT 
因 测 序 和 质谱 分 析 进 


WEA) 18 个 差异 表达 基 
， 但 已 经 确定 是 病程 相关 和 蛋 
EAN a 14 THS 的 3 个 差异 基因 ( Cluster-24362.0, 
侵 染 后 其 表达 量 都 会 升 高 ， 且 在 qRT-PCR 验 
关 ， 后 期 可 着 重 分 析 THS 酶 的 


酶 与 抗 根 肿 病 密切 相 


AAR 


HA FSG MOA 
PDE EME AI RAK AE 


表达 量 上 调 ， 这 是 因 
匆 对 北 境 的 耐 受 性 
因 组 ， 无 法 对 中 
BOLE 


Al 


有 


步 探索 完整 的 叫 喉 生物 碱 合成 通 


KITE EME RA 


PR-10 超 家 族 成 员 之 


THRE BRA RR EFA YA HT AS, ECHR SPR Ee SE (TAT) 能 催化 酪 氨 酸 合 
成 4- 羟基 苯 丙 酮 酸 ， 现 已 在 拟 南 芥 (Lopukhina etal., 2001) . F2 (Huang et al., 2008) 、 
紫 苏 《上 昌 晓 玲 等 ，2012) 等 植物 成 功 分 离 克 隆 了 TAT 基因 ， 并 发 现 植物 激素 水 杨 酸 和 脱落 
酸 处 理 后 可 在 转录 水 平 上 诱导 该 基因 的 表达 。 本 研究 结果 也 发 现 编码 TAT 的 基因 
Cluster-31730.0 和 Cluster-28040.6 在 根 肿 菌 胁迫 后 表达 量 上 升 ， 说 明 TAT 酶 在 胁迫 中 发 挥 
重要 作用 。 

qRT-PCR 实验 结果 表明 , 在 根 肿 菌 胁迫 下 ，3 个 时 间 点 称 蓝 实验 组 与 对 照 组 相 比 ，7TZEy 
(BER, 2020) ~ AUX/IAA (Jie et al., 2018) ~ GH3 (Lu et al., 2022) 、S4UR (Li et al., 2022). 
PYR/PYL (Kim et al., 2020) ~ PP2C (Yu etal., 2019) 、ABF〔 陈 乃 钰 等 ，2021) . B-ARR 
(Falconieri et al., 2022) ~ TGA (Qi etal., 2022) ~ AHK (Cerbantez-Bueno et al., 2020) 这 些 
基因 表达 量 都 上 调 且 随 着 根 肿 菌 侵入 时 间 的 延长 , 同一 基因 其 表达 量 不 断 降 低 , 准确 反映 了 
根 肿 菌 胁迫 下 大 多 数 基因 在 徐 蓝 中 的 表达 模式 。 

综 上 所 述 , 对 根 肿 菌 侵 染 的 获 蓝 进行 转录 组 测序 和 代谢 组 分 析 , 极 大 的 丰富 了 根 肿 菌 胁 
迫 下 其 蓝 的 生物 学 信息 , I Be EG] RSS AE YO Ti AL SFE IE RS AE YO TBS CE BE A], 探讨 
了 这 些 关 键 基因 在 应 对 逆境 胁迫 下 的 表达 规律 , 为 后 续 深入 研究 这 些 基因 的 功能 、 解析 根 肿 
菌 胁迫 下 获 蓝 生物 碱 的 积累 机 制 黄 定 基础。 
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